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Abstract 

We describe key elements of the perturbative similarity renormalization group procedure for 
Hamiltonians using two, third-order examples: 0^ interaction term in the Hamiltonian of scalar 
field theory in 6 dimensions and triple-gluon vertex counterterm in the Hamiltonian of QCD in 4 
dimensions. These examples provide insight into asymptotic freedom in Hamiltonian approach 
to quantum field theory. The renormalization group procedure also suggests how one may 
obtain ultraviolet-finite effective Schrodinger equations that correspond to the asymptotically 
free theories, including transition from quark and gluon to hadronic degrees of freedom in case 
of strong interactions. The dynamics is invariant under boosts and allows simultaneous analysis 
of bound state structure in the rest and infinite momentum frames. 
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1. INTRODUCTION 



Canonical Hamiltonians of quantum field theories can be written in terms of creation and 
annihilation operators. Let us denote those operators by q. In this notation, for example, a 
triple-gluon vertex in the light-front QCD Hamiltonian has the structure 

Hy = ^ f[123]6il + 2-3)gY,23qlqlq3 + h.c. , (1.1) 

123 

where symbols 1, 2 and 3 denote colors, spins or momenta of gluons, [123] is a shorthand 
notation for the integration measure over the gluon thrcc-momcnta, S{1 + 2 — 3) is the Dirac 
(5-function of three-momentum conservation, and Y123 is a function of the gluon quantum num- 
bers, implied by QCD. g denotes a bare canonical coupling constant, which is expected from 
Lagrangian approach to require ultraviolet renormalization. However, the Hamiltonian lacks 
many of the Lagrangian density symmetries that are employed in the perturbative Lagrangian 
renormalization procedure and, in all light-front Hamiltonians, transverse and longitudinal di- 
rections are treated differently. For that reason, and to introduce a method for renormalization 
of Hamiltonians such as (1.1), we need to review the origin of the ultraviolet divergences. 

One can evaluate matrix elements of Hy between states of the form \ijk...) = g|gjg^...|0). 
Consider (12|iJy|3), which is proportional to Y123. The trouble is that Y123 does not vanish 
when the relative motion of particles 1 and 2 becomes very energetic. In other words, the 
interaction Hamiltonian directly couples states of small kinetic energy to states of arbitrarily 
high kinetic energy. For example, when we square Hy and evaluate I]j -f^y |«) («|i^y, the range 
of \i)s to sum over, on the energy scale, is infinite. The sum diverges, and the square of Hy 
does not exist. If one tries to find eigenstates of the Hamiltonian, the divergence will dominate 
finite terms. Also, exp —iHt diverges and no conclusions can be drawn from knowing H as it 
stands. The problem is worse than the inifinite energy range of interaction implies by itself: 
the function Y'123 in QCD grows when the energy difference between kinetic energies of particles 
1 and 2, and particle 3 grows. The larger are the kinetic energies of the intermediate particles, 
the more important become the interactions, and we are sent into an abyss of Fock space states 
without bounds. All physically relevant local quantum field theories have this trouble. 

To see the essence of the problem, imagine a matrix of the Hamiltonian matrix elements in 
the basis of eigenstates of certain Hq, HQ\i) = EQi\i), H^j = {i\H\j). Our problem is that the 
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corners marked low — high and high — low on the above figure, contain too large matrix 
elements for the Hamiltonian matrix to have eigenvalues that are independent of the matrix 
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boundaries. In order to understand the boundary-dependence problem of the eigenvalues, we 
put an upper bound, denoted by A, on the basis states energies and we work out what happens 
with eigenstates of Hamiltonian matrices whose size is limited by the conditions E^i < A and 
Eqj < A, when we change A. In particular, we ask what properties must Hy have for the 
spectrum of H/\ to have a limit when A ^ oo. 

Wilson asked this question in case of a model Hamiltonian with big energy gaps between 
successive energy scales that were included in his calculation, and he studied the influence 
of coupling between small and large energy states on the lowest eigenvalues. His method 
for dealing with the A-dependence of the spectrum in the multiscale eigenvalue problem (in 
general, the problem of dependence on regularization, of any kind), was based on an iterative 
procedure. Initially, one solves the highest energy part of dynamics and focuses on its lowest 
eigenvalue levels whose dynamics, in turn, is dominated by states with energies lower by one 
energy gap. Then, one solves this next lower energy scale dynamical problem, and one repeats 
the process many times. Starting from the energy scale A, one eventually arrives at a finite 
scale A. This is schematically indicated on the following figure. 

Wilson's approach 



A 



Rx 




A 



In this figure, we see a new small matrix of size A, denoted by Hx- This matrix is calculated 
using an operation Rx, which is constructed in the sequence of steps lowering the cutoff from 
A to A. The construction is based on the principle that the smallest eigenvalues of the small 
matrix, Hx, should be the same as the smallest eigenvalues of the big matrix H^. The algebraic 
derivation of the small Hamiltonian Hx is designed to guarantee the equality of the smallest 
eigenvalues. The crux is that if all matrix elements of Hx are independent of A (in general, 
independent of the regularization one uses to define Ha) then, the eigenvalues of Hx must be 
independent of the regularization. Therefore, if we know what to do with the regularization 
dependence of matrix elements of Hx, then we know how to go about regularization dependence 
of the spectrum of H^- 

Note, that A is finite and can be chosen arbitrarily, as long as A ^ oo. The set of transfor- 
mations that connect Hamiltonians with different values of A is called renormalization group. 
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Physical results should be independent of A, by construction. It is clear that H\ cannot be 
equal to merely PxH^Px, where Pa denotes a projection operator that projects on the space 
of states with energies Eq^ < A. Some additional terms must be included, which reproduce 
dynamical effects from above A. Similarly, we do not expect that is merely PAHcanonicaiPA, 
where H canonical is built from Hq and terms such as (1.1). By the same argument as for Hx, 
some additional terms are required in Hj\, to include effects from above the cutoff A. We 
will call those terms counterterms. The problem is how to find them. According to Wilson, 
they are found from the condition that all matrix elements of Hx become independent of the 
regularization in the limit A — oo, for all finite values of A. One must take care not only of 
the diverging (i.e. A-dependent) regularization dependence, but also of the finite regularization 
effects. The coupling constant g in Eq. (1.1) is changed in the renormalization process as a 
result of introducing counterterms. 

The problem with the transformation Rx is that the large energy gaps are absent in phys- 
ically relevant cases, and perturbation theory based on energy scales alone fails. One can try 
to take advantage of a small coupling constant but we know that naive perturbative expansion 
does not work in degenerate cases. For example, one may think of effects familiar from elemen- 
tary degenerate perturbation theory for eigenvalues of a Hamiltonian matrix with a few rows 
and columns. We know that perturbation theory cannot work unless one properly chooses the 
initial basis states in the degenerate subspace - eigenstates of the interaction matrix (if it is 
a few rows and columns). Only in that basis the perturbative limit g ^ exists. Otherwise 
perturbation theory produces vanishing energy denominators that lead to diverging terms and 
the calculation is misleading. In case of quantum field theories of interest, the situation is much 
more involved than in simple matrix case due to multiply degenerated continuous spectra of 
Hq. Moreover, in asymptotically free theories, we expect that the interaction strength grows 
when we go from H/^ to Hx and small energy denominators are certainly expected to produce 
large effects. Thus, the degeneracy of spectra and strength of couplings do not allow us to do a 
precise analysis of Rx and Hx- So, the operation Rx is of limited applicability in the canonical 
approach to quantum field theory. 

2. SIMILARITY FOR HAMILTONIANS 

There is an alternative approach (see the figure below). Instead of calculating a 

small Hamiltonian matrix, we can also calculate a narrow matrix. Namely, a similar (in the 
sense of algebraic similarity) matrix that has the same eigenvalues but whose matrix elements 
Hxij vanish if \EQi — EqjI > A (or another condition of "narrowness" is satisfied - hermitian 
matrices can be diagonalized and, therefore, partial diagonalization to a narrow matrix should 
be possible). The choice of near-diagonal form is motivated by the following property of near- 
diagonal matrices: when we act with them on a state of some finite energy, a single action of the 
matrix can rise the energy by at most A, i.e. by its width on the energy scale. In perturbation 
theory for eigenvalues of a near-diagonal matrix, corrections will not be sensitive to the cutoff 
A up to the order n ~ A/(2A), since one has to go up in energies and come back to the initial 
energy range through action of interactions, and at least n are needed to reach the boundary 
starting from A. In similarity, the crux is that if matrix elements of the narrow matrix of finite 
width A are independent of regularization when A — > oo, then the spectrum of Hx will be 
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independent of regularization to all orders of perturbation theory. 

Similarity approach 




The transformation Sx is called similarity transformation. The effective Hamiltonian ma- 
trices with various widths A are connected by transformations that are called similarity renor- 
malization group transformations for Hamiltonians. We limit our discussion to unitary trans- 
formations S\. The key feature of the similarity approach is that perturbative construction of 
S\ avoids small energy denominators entirely - they are limited from below by the width A. 
In turn, the perturbatively calculated narrow Hamiltonians can be diagonalized numerically, 
which is the ultimate way to find solutions to complex non-perturbative problems of the original 
theory. One could ask, why don't we go all the way to A = 0, which would mean complete 
diagonalization through S\ with A = 0? This is impossible in perturbation theory. We can 
trust perturbation theory for calculating Hx only for not too small values of A. 

In other words, the perturbative similarity transformation S\ involves energy changes that 
are at least as large as A and the problem with large effects in perturbative evaluation of effective 
Hamiltonians is overcome. But that does not mean we eliminated any of the nonperturbative 
effects. They are still hidden in the narrow effective Hamiltonian, as much as they were in 
the initial one. The only thing we accomplish through similarity, is the elimination of direct 
couplings between states of interest to us and very high energy states. This is a prerequisite 
that we need to define an ultraviolet finite, nonperturbative Hamiltonian eigenvalue problem 
in quantum field theory. 

One can apply various perturbative procedures for calculating Sx and Hx- Wegner invented 
a beautifully simple scheme for evaluating near diagonal Hamiltonians in solid state physics. 
It was shown that Wegner's equation can be employed in the renormalization group scheme. 
A number of variations reported in the literature in different areas, is growing. 

The bottom line is that when solving for the spectrum of a narrow Hamiltonian, we do not 
have to diagonalize the whole matrix of size A ^ oo. We can select a window, as is illustrated 
in the next figure. Diagonalization of that small window is much simpler than diagonalization of 
the whole matrix and the window eigenvalues match the whole matrix eigenvalues in the middle 
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range of window energies. [§] The reason is that the wave functions have width comparable 
to A on the energy scale, which is indicated on the figure (see also Appendix B in Ref. 0). 
One cannot expect wave functions of eigenstates of the initial Hamiltonian matrix, i^Aj to 
be dominated by some small range of energies. In contrast, the eigenstates of matrix Hx are 
expected to have wave functions that have this property. Recent calculations of quarkonium 
and glueball spectra in hght-front QCD exploit this feature. 0] [|T^ 



Window Hamiltonians 



Sy 




Ref. P outlines the similarity procedure starting from an initial Hamiltonian, through 
evaluation of the narrow Hamiltonians (having found the necessary counterterms in the initial 
Hamiltonian) in perturbation theory, to diagonalization of a small window to get the bound 
state energy, using a matrix example. The matrix model is asymptotically free and has a 
bound state. The coupling constant is a function of the effective Hamiltonian width A, we say 
it "runs". For example, it may equal about 0.06 at A = 65 TeV and about 1 at 1 GeV. Still, the 
window Hamiltonian of a few GeV size can be calculated in second order perturbation theory 
and the window bound state eigenvalue deviates from exact solution by only 10%. Once we 
understand that example, we can return to quantum field theories with interactions of the form 
(1.1) and attempt a calculation of the corresponding "window" Hamiltonians, with running 
couplings. Calculations can be carried out using the notion of effective particles. 



3. SIMILARITY FOR PARTICLES 

When we deal with huge Hamiltonian matrices of quantum field theory the number of states 
is as big as we let it be and the number of matrix elements we have to think about becomes very 
quickly incomprehensible. We have to reduce the amount of information that we need to know 
at the beginning. Imagine we would know matrix elements of the interaction a/r in atomic 
basis functions, numerically, but we would not know that they all correspond to the Coulomb 
force. It would be very hard to relate what happens in one atomic system to what happens 
in another one. Therefore, when we aim at universal calculations of effective Hamiltonians in 



6 



theories that contain interactions such as (1.1), we may proceed to a new version of similarity 
transformation, which avoids deahng directly with Hamiltonian matrix elements in a particular 



basis and, instead, operates at the level of field operators. |jTT]] 

Let us introduce a unitary transformation U\ that transforms field operators (denoted here 
by 0, independently of their spin or other quantum numbers they carry), 

(t)x{x) = Ux(t)oo{x)Ul . (3.1) 

(j)oo{x) denotes a bare quantum field operator that, at any prescribed time, can be expanded 
into creation and annihilation operators for bare particles in a canonical fashion that we do 
not need to define here very precisely. (f)x{x) denotes an operator that is built in exactly the 
same way from creation and annihilation operators for effective (dressed) particles. This kind 
of transformation is motivated by physics of hadrons, whose structure can be explained in 
a constituent quark model. Dressed particles in a given theory interact differently than the 
bare ones. Namely, bare ones have interactions like (1.1), while the effective ones can only 
exchange momentum transfers that are limited by A. This is secured by the construction of 
Ux, to be explained below. Therefore, the effective particle wave functions of eigenstates of 
the Hamiltonian may quickly fall off when momenta or number of the effective particles deviate 
from the physically dominant values. This is why one can hope to obtain a constituent picture 
of hadrons in QCD using similarity for particles. More generally, the expected convergence of 
eigenstate expansion in effective particle basis in Fock space opens a door to studies of few-body 
systems in quantum field theory. 

In order to set up equations that will allow us to calculate Hamiltonians for effective parti- 
cles, let's rewrite Eq. (3.1) in terms of the creation and annihilation operators, 

qx = UxqooUl . (3.2) 

All we need to do next is: take the bare Hamiltonian of our theory, as it is given initially in 
terms of goo, calculate counterterms it needs to contain in addition to the canonical terms, 
obtain this way our initial -ffoo(9oo), and rewrite it in terms of q\. U\ is secured to be unitary 
by construction. The whole point of the construction is that the resulting H\{qx) is to contain 
only such interaction terms that, when we evaluate their matrix elements between Fock basis 
states of effective particles, the resulting effective Hamiltonian matrix is narrow, of width A, as 
in the similarity procedure for Hamiltonian matrices we discussed in previous Section. It will 
not be necessary to go into details here. Only a brief outline of the scheme follows. 

Since rewriting the Hamiltonian in different degrees of freedom does not change the operator 
itself, we have H\{q\) = Hoo{qoo)- One may think about H\{q\) as a QCD Hamiltonian written 
in terms of constituent quarks and gluons, and about i?oo(goo) as the same QCD Hamiltonian 
written in terms of canonical quarks and gluons, associated with partons, or current quarks 
(to make the connection between hadronic rest frame constituents and partons in the infinite 
momentum frame, we have to use the light-front form of Hamiltonian dynamics, see |T^ for an 
outline of light-front QCD in the context of renormalization group procedure for Hamiltonians). 

Applying the transformation U\, one obtains Hx = Hx{qoo) = U\Hoo{qoo)U\- This relation 
means that the operator 'htx has the same coefficient functions in front of products of q^o as the 
effective Hamiltonian Hx has in front of the unitarily equivalent products of qx- Differentiating 
Tlx one obtains 
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A^Hx = - [r,,n,] , (3.3) 

where the generator Tx is related to Ux by 

= . (3.4) 

The script letters are introduced to indicate that the operators can be conveniently thought 
about as expanded into sums of products of operators goo- The latter are independent of A and 
are not differentiated in Eqs. (3.3) and (3.4). In other words, Eqs. (3.3) and (3.4) describe only 
the flow of coefficients in front of the creation and annihilation operators. Effective Hamiltonians 
are obtained from TCx using Hx{qx) = UxHxU'l- 

The key element now is how one defines Tx. This is the domain of similarity for effective 
particles. In its essence [Q, one studies what one has to do to get the narrow Hamiltonian 
matrices as a result of the procedure, and these studies tell us what to put for Tx- There exist 
infinitely many choices. The one that the present author used to get results described in the 
next Sections is of the following form [11| [1^ 



[Tx,n,x] = ^(1-^a)[^a] . (3.5) 

The symbols Q and F require explanation. The effective Hamiltonian Hx contains form factors 
of width A in all its vertices. If we denote an operator without the form factors by Gx, our 
Hamiltonian takes the form Hx = Fx[Gx\, where the operator Fx inserts the form factors. 
With these form factors, momentum transfers in interactions between effective particles are 
guaranteed to be at most of the order of A. Qx = l^lfixl^x- We divide Qx into two parts, a part 
that is bilinear in q^, and an interaction part that would vanish if the coupling constant were 
equal 0, so that Qx = Qq + Qix- The operator Qjx satisfies the following differential equation as 
a consequence of Eqs. (3.2)-(3.5), 



/^.Ud-/)^; 



(3.6) 



We dropped the subscript A on the right-hand side for clarity. / denotes the similarity form 
factor introduced by Fx and the curly bracket with the subscript Qq denotes a solution for Tx 
resulting from Eq. (3.5). 



4. ASYMPTOTIC FREEDOM IN SCALAR THEORY 

Since the interaction term (1.1) is only a part of the QCD Hamiltonian and the function 
Yi23 depends on spins and momenta of gluons, let us first discuss the case of scalar field with 
classical Lagrangian density 

/: = ^(9^09'^0-/iV)-|0' • (4.1) 

In this case, the interaction term in the corresponding Hamiltonian is of the form (1.1), but 
^123 = 1/2 and calculations are much simpler than in QCD. Our goal is to describe results for 
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the light-front Hamiltonian for effective bosons calculated in perturbation theory up to third 
power in g. Although our presentation is based on Ref. |jl3[ that uses plain expansion in powers 
of the reader may also wish to compare our results with Ref. [Q, where a different scheme 
is used, including transverse locality and coupling coherence. |T0[ 

The light-front Hamiltonian corresponding to the Lagrangian density (4.1) reads 

/fc"*" ^ -|- ^ 
[k] a^fc«oofc + 

(4.2) 

where is a smooth regularization factor and Xa denotes counterterms (derivable in pertur- 
bation theory). In n dimensions, [A;] means 6{k^)dk~^d'^~'^k-^/{2k^{2'n-)"'^^). We choose 



= exp -''"+J''^" , (4.3) 



where xi = kf /k^ and Kj'-g = k^ — xik^, rji = rj{xi), and ?7 is a useful function of its argument. 
A natural choice is ri{x) = 1, for it is simple. Leaving r] unspecified will help us identify finite 
regularization effects. 

The similarity form factor for an operator containing u creation operators and v annihilation 
operators is defined by 

Mu,v) = exp[-iMl-Mlf/\'] . (4.4) 

The script notation for invariant masses means Ai^^ = {ki + ... + A;„)^, where the minus com- 
ponents of the momentum four-vectors are given by k^^ = {k^"^ + iJ.'^)/kf for i = 1, and 
similarly for v. 

Equation (3.6) can now be solved order by order using expansion in powers of g. Firstly, 
one obtains the counterterms Xa as the initial conditions at A = cxd that render regularization 
independent finite A Hamiltonians. To order g^, the regularization dependence of Hx lets us 
identify two counterterms: the mass counterterm 

,2 



/3ooii = / [k] -^aLfcfloofc , (4.5) 



and the vertex counterterm 



7oo2i = [kik2k3]2{2TT)^S^{ki + k2- ks)'yoo{ki,k2,k3)a^^^y^,^^aook3rA ■ (4.6) 



Without loss of generality, we assume that some gedanken experimental data require the 
mass squared parameter in effective Hamiltonian with A = Aq to be equal /i^ + Sfi^. This means 
that when one calculates observables using the effective Hamiltonian, must equal fi'^ + 6nl 
to fit the data. This condition, by tracing the renormalization group equation for Hx back to 
X = oo, tells us that 
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2 J 2(27r)5 Jo x{l-x 



dx 



rA0 



(4.7) 



The script Ai denotes invariant mass, Ai^ = (k^^ + /i^)/x(l — a;), and the regularization factor 
is 



Tap = exp[-2[r]{x) +7]{l-x)]K^yA'^} 



(4.8) 



Integration gives two diverging terms, one proportional to and another one proportional 
to log A. The remaining finite part depends on our choice of the function rj. For example, 
evaluating the integral for ri{x) = 1/x one obtains 



6 fx. 



2 

oo 



(4vr) 



24 



rA^ 



25 , A 2 
^ - log — + 
6 /i ' 



(4.9) 



where /i^ has a finite limit when A ^ oo. The logarithmically divergent part is independent of 
the function r] and agrees with results for the Lagrangian mass squared counterterm obtained 
using Feynman diagrams and dimensional regularization |jl5| |jT6| in the following sense: when 
one changes A to A' the logarithmic part of the counterterm changes with A as the mass 



squared changes as a function of the renormalization scale in Eq. (7.1.22) in |jT6|. 

The vertex counterterm is defined by the requirement that the effective vertex in the Hamil- 
tonian is free from regularization dependence for arbitrary finite values of Aq. The one loop 
regularization sensitive contributions to the effective vertex function are given by 



X 



2(27r)^ 



X 



dx 



1 x(l — x){x — Xi) Jo 
1 dx 



2 J 2 o ^1 

K dn 8- 



x[l — x) Jo 



oo _Q 
2 J 2 

K dn 



exp 



XX2J^^ 



exp 



"A2~ 



X2) 



+ 



(4.10) 



where 



Cr, = r]{x) +-r]{l-x) + {r]{xi/x) + ■r][{x - Xi)/x\} {xi/ x)"^ + r][{x - Xi) / X2]+'r][{l - x) / X2] (4.11) 
and 

dr, = 2[r]{x)+r]{l-x)] . (4.12) 

The counterterm function 700(^1,^2,^3) must remove the regularization dependence from the 
above expression. The regularization effects are independent of Ki2- Dropping all parts that 
are independent of regularization, we conclude that 
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1 

(47r)3 



A 

3 log 4 



1 1 — X 
(ix logc^ + [xi ^ X2) 

XI X2 



+ 3/ dx x{l — x) log dr. 
Jo 



(4.12) 



must be independent of regularization. We see that the diverging regularization dependence 
of the interaction vertex, i.e. the term proportional to log A, is independent of the particle 
momenta and one can remove the divergence by introducing a 7oo(^i, ^2, ^3) that is equivalent 
to changing the initial coupling constant g in Eq. (4.2). Thus, no diverging x-dependent 
counterterms are required - a different situation than in |l^. However, it is also visible that 
the vertex contains a finite regularization dependent part that is a function of Xi. The function 
depends on our choice for 77. For example, if 77 = 1 one has = 4 + 2(xi/x)^ and drj = 4. The 
resulting integral is a function of xi, and needs to be subtracted. But this would not assure 
us that the whole ultraviolet regularization dependence is removed, because we work with a 
specific functional form of the regulating function (4.3). 

Since the whole regularization effect is independent of A and it can be completely 
removed from the effective interaction by subtracting its value for = at an arbitrarily 
chosen finite Aq. However, one has to add back the finite regularization independent part of 
the effective vertex, which is a function of xi, denoted below by 70 (xi). The function 70 (xi) 
is necessary to recover Poincare symmetry of observables, because our regularization spoiled 
the symmetry. The symmetry may be restored once counterterms remove the regularization 
effects, but one is not allowed to change terms independent of the regularization, which were 
given by the initial covariant Lagrangian density unambiguously. Therefore, the function 70(2:1) 
must be reinserted. This function is not altered when A changes and could be considered 
marginal in analogy with usual renormalization group analysis. The ultimate adjustment of 
the function 7o(xi) requires 4th order calculations. For there exists in 0^ theory no 3rd order 
scattering amplitude one could use to find out what function 70(2^1) renders Poincare symmetry 
of scattering observables with our choice of ta in Eq. (4.2). However, it should be pointed 
out that the function does not influence the way the 3rd order running coupling constant in 
effective Hamiltonians depends on A. 

Thus, in Eq. (4.6), the counterterm function 700(^1,^2,^3) = 7oo(a^i, ^12), "which removes 
the regularization dependence from the effective vertex reads 



7oo(xi,K]^2) = -7Ao(a;i,0^) + 7o(a;i) • (4.13) 

This result is used to define the new regularization dependent coupling constant g/\ in the initial 
Hamiltonian in Eq. (4.2). We select a convenient value of Xi = Xq and obtain 



9a 9,1 n-L^ 



I - 7Ao(a;o,0^) +70 



(4.14) 



where 70 = 7o(xo). We see that the initial coupling g is replaced by a new A-dependent quantity 
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+ o(/) , (4.15) 

with certain free constant niQ. Thus, the theory exhibits asymptotic freedom in 3rd order 
terms. Our result agrees with hterature, say Eq. (7.1.26) from []rB|, in the sense that when we 
change A, the change required in the couphng constant in the initial Hamiltonian for obtaining 
A-independent effective Hamiltonians matches the change implied by Feynman diagrams and 
dimensional regularization. 

Having established the structure of counterterms we can proceed to evaluation of the finite 
similarity flow of effective Hamiltonians towards small widths A. The effective kinetic energy 
term in narrow Hamiltonians is 

7 _L 2 2 

Hxu = J [k] -^^ala,, , (4.16) 

where 



9A = 9 



1-9' 



4(47r)' 



log- 



A 

mo 



2 2 I r 2 

(4.17) 

The above result is particularly simple for /i = and in that case it reads {Sfi^ is proportional 
to g^) 

Logarithmic dependence on A arises for /i > 0. The value of 6fil could be found, for example, 
by solving a single physical boson eigenvalue problem, expressing the physical boson mass in 
terms of 5/iQ and adjusting the latter to obtain the gedanken experimental mass value for bosons. 
Note a change in the mass function of cutoff parameter, from the case of the mass counterterm, 
dependent on A, to the case of running mass term, dependent on the width A (independent of 
A). The change corresponds to a transition from the initial side of a fixed point (bare canonical 
Hamiltonian with regularization) to the other side of the fixed point (renormalization group 
trajectory of effective Hamiltonians in the similarity procedure, cf. 0). 
The effective vertex reads 



Hx2i = J [kik2h] 2{2T:f5\k^ + k2 - k-^) fx[{ki + k2)\ kl] Vx{xi, k^^) a{^^a{^^axk, , (4.19) 

where Vx{ the effective vertex function and fx is the similarity vertex form factor. 

The vertex function is given by an integral over loop variables x and k-^ of a known function. 

m 



12 



We define the running coupling constant as the value of 2V\{xi, /tfg) at a chosen configuration 
of momentum variables, denoted by {xiq, k^2o)- other words, gx = 2Va(xio, i^uo)- possible 
choice for massless bosons is Xio = and K^^go ~ 0- This is a natural definition, analogous to the 
standard Thomson limit in the case of electron charge in QED. This choice greatly simplifies 
the integrand, giving Va(0, O"*-), so that the result can be fully produced here (5^0 is the value of 
gxo required by phenomenology done using -f^Ao) 



3 1 1 r°° dz 
9^ = 90 + 9o^J^-sl ~ 



2(A - fi) - 2(/o - m + 2QUi - fx) - nfo - fo) + 9(/o^ - /, 



(4.20) 

where f\ = exp —z^jX^ and /o = exp —z^jX^. A straightforward integration gives 

9. = 9.- ^0^(1^ log ^ ' (4-21) 

which exhibits asymptotic freedom. Differentiating with respect to A and keeping terms up to 
order g\ one obtains 

This equation demonstrates the same /? function for coupling constants in effective Hamiltonians 
as obtained in Lagrangian approaches using Feynman diagrams and dimensional regularization, 
when one identifies the renormalization scale with the Hamiltonian width A. This is encour- 
aging but one needs to remember that for comparison of perturbative scattering amplitudes 
in Hamiltonian and Lagrangian approaches it is necessary to make additional calculations and 
at least of fourth order in g. Beyond model matrix studies such as in 4th order similar- 
ity calculations have so far been carried out only in a simplified Yukawa model by Maslowski 
and Wi§ckowski [|l^ (the latter model calculations should be helpful in setting up a light-front 
theory of nucleons and pions). 

Integrating Eq. (4.10), one obtains [a = 5'^/47r) 

l + ao(3/327r2)logA/Ao ' ^^'^^^ 

which shows our result for a boost invariant running coupling constant in effective Hamiltoni- 
ans. Our procedure explains how the running coupling constant can be included in quantum 
mechanics of effective particles, which is given by the Schrodinger equation with the corre- 
sponding Hamiltonian Hx. One can evaluate matrix elements of the small width Hamiltonian 
in a limited subspace of Fock states built of the effective particles. The form factor in the 
interaction vertex (4.19) secures a small range of the interactions on the energy scale and one 
can expect a rapid convergence of wave functions in the effective particle basis. 

Note three characteristic features of the Hamiltonian calculation. (1) No field renormaliza- 
tion constant appeared, since the similarity transformation did not eliminate (or integrated out) 
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any degrees of freedom. (2) No vacuum effect played any role, since we used the light-front form 
of dynamics. Extensive literature concerning the vacuum issue can be traced through reference 
p!9| . 0^ theory is unstable due to a possibility that the field takes an infinitely large, negative 



value. It would be interesting to check if the perturbatively evaluated effective Hamiltonians of 
small widths have any tendency to develop eigenstates that deviate in that direction. (3) The 
Hamiltonian structure is invariant with respect to boosts, including boosts from the rest frame 
of any bound state to the infinite momentum frame. This suggests the approach outlined above 
should be tried in QCD, and in effective theories of strong interactions in nuclear physics, to 
connect low energy observables, such as binding energies, radii or magnetic moments of bound 
states, with high energy ones, such as parton distributions, form factors or jets. 



5. QCD GLUON VERTEX COUNTERTERM 

We come back to Eq. (1.1) in QCD and repeat the same analysis as we did for the scalar 
theory. Most of the procedure remains the same, but an important complication arises. The 
vertex function in the canonical Hamiltonian has now the form (c refers to color and e to 
polarization of gluons) 

1^123 = tr^'^n^^e;^ ■ ei^i, - el\i ■ €^^^2/^2 - e^e^ ■ el^n^Jx,] , (5.1) 

in which the characteristic factors of /x tend to infinity when x — >■ 0. Ultraviolet coupling 
constant divergences in 0^ theory in 6 dimensions resulted from transverse momentum inte- 
gration / d^K^/K'^, where I/k^ came from the two denominators of third order perturbation 
theory. In QCD in 4 dimensions, we have instead / (pK-^{n-^ / xY / k^. Therefore, in QCD (more 
generally, in gauge theories), we have to introduce a separate regularization of small x behavior 
of interaction vertices in the Hamiltonians. For example, in the QCD counterpart of Eq. (4.2), 
we have to insert a factor, denoted by r^, in addition to ta, that will effectively cut-off the 
region where one of the gluons 1 or 2 carries a smaller fraction of ^3" than the size of a small 
parameter 5. 

The small x regularization function may appear to be only a technical detail. But it 
was pointed out by Perry that singularities at small may be related to effective confining 



potentials in Hx^ calculable already in second order perturbation theory. In short, the 
canonical light-front QCD Hamiltonian contains terms that are singular at small and the 
singularity contributes to the effective Hamiltonians H\, providing potentials that grow with 
distance between color charges. This is quite different a situation from other formulations 
of the theory, where second order calculations are not expected to tell us anything about 
confinement. Therefore, the small x features of QCD in the light-front Hamiltonian approach 
deserve extensive studies. Here, we merely report some initial results for third order gluon 
vertex counterterm, indicating x-dependent features. 

The whole analysis of the previous Section can be repeated step by step and one can derive 
the interaction term for effective gluons, in the narrow Hamiltonian Hx for QCD. The condition 
that the effective vertex is independent of regularization parameter A gives us the diverging 
triple-gluon vertex counterterm in the initial QCD Hamiltonian. The new features appear in 
the coefficient of log A. Namely, the diverging part of the vertex, to be compensated by the 
ultraviolet counterterm, has the form 
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£^ [^iVc-^% + iVJ(a;i,X2)] logA-Yi23 , (5.2) 

where the function /(xi, X2) is symmetric in its arguments. This function originates from three 
successive actions of the triple gluon interaction, the gluon mass correction providing only a 
constant contribution, and it depends on the regularization factor in the initial Hamiltonian. 
In fact, /(xi,X2) = [-logxi+/o rfx(2/x+l/(l-a;))(r53-r52) + (l ^ 2)]/2, where r^g = rs3{x,Xi) 
is a product of three factors rs for the three successive triple-gluon interactions, and rs2 = rj{x), 
is a product of two factors rs from the gluon mass counterterm. For rs{x) = 9{x — S)9{1 — x — 6) 
in the bare vertex, for gluons carrying x and 1— a; of the single gluon momentum in Eq. (1.1), 
f{xi,X2) = log [mm (xi, 0:2)], which is negative and hence reduces the rate at which the initial 
coupling depends on A. The latter feature comes about as follows. Since the counterterm must 
contain the term opposite in sign to (5.2), the coupling constant in the regularized ifoo(Q'oo) in 
QCD is changed to [cf. Eq. (4.15)] 



9A = 9 - 



9^ rll 1 

—Nc - -Uf + NJ{xo, 1 - xo] 
. 1/ 



log A . (5.3) 



47r2 

The coefficient of log A is independent of the parameter 6, but it depends on the small x 
regularization in a finite way. Having derived the function f{xi,X2), for certain choices of 
Xi = Xo in the definition of the Hamiltonian coupling constant, one could obtain triviality 
instead of asymptotic freedom. This is a peculiar result for the regularization given in Eq. 
(4.3) and the sharp cutoff on x at 6. In the limit when both rg^ and rs2 are replaced by 1, 
one would obtain /(xi,X2) = — log ^/xiX2 , which is always positive, and would accelerate the 
asymptotic freedom rate of change of g/\ with A. One can seek choices of regulating functions 
ta and rs that eliminate the unusual logarithm (the other terms are standard) but it is not 
known if a finite function of Xi is not necessary in place of /(xi, X2) in light-front Hamiltonians 
anyway, to restore symmetries for physical quantities. In addition, as in the scalar theory, the 
ultraviolet finite part of the counterterm involves an unknown function of Xi. One could say 



that the structure of the model from Ref. [1/] is closer to QCD than to scalar theory in 6 



dimensions. Evaluation of the effective coupling constant gx in QCD, in analogy to Eq. (4.20) 
in scalar theory, may shed some light on how to disentangle genuine ultraviolet from small x 
singularities in Hamiltonians. 

It is clear from the above example that effective light-front Hamiltonians of QCD require 
careful studies employing various types of regulators before we will know the optimal ways 
of calculating window Hamiltonians. The interplay of transverse and longitudinal momentum 
variables may lead to surprising results. However, the calculations are certainly doable and 
the resulting matrices will tell us about details of QCD dynamics in the Fock space of effective 
quarks and gluons. The similarity renormalization group procedure for Hamiltonians is able 
to reveal new features of effective particle dynamics which standard Lagrangian approaches do 
not reveal. 

6. TRANSITION TO NEW DEGREES OF FREEDOM 

We may hope to make a transition from the effective QCD degrees of freedom to nuclear 
physics hadronic interactions, such as pion-nucleon coupling, after we achieve understanding of 
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the narrow Hx in QCD. To see the basis for this hope, let us come back again to the matrix 
picture. 

Once we have the narrow Hamihonian matrix, we can divide it into the boxes as it is 
illustrated in the figure below, neglecting the small triangles in the band outside of them, 
which are initially left out. We can find eigenstates of the boxes (they correspond to different 
invariant mass states) and calculate the band diagonal matrix matrix elements in the basis built 
from those eigenstates. The states corresponding to the middle energy scale of each box will 
not interact very strongly with neighboring (in energy, or mass) states, but matrix elements 
sensitive to the left-out triangles will lead to strong interactions. We can imagine that the lowest 
box corresponds to nucleons interacting through potential forces, the next box corresponds to 
nucleons plus one meson, the second box to nucleons and two mesons, etc. This is how one 
could make a connection between the QCD dynamics and nuclear physics through similarity 
renormalization group for Hamiltonians (cf. ||^1|). 



Changing degrees of freedom 




However, this is not the only possibility one can try to explore for the change of basis. 
One can consider new basis states built of quarks and gluons, possibly open gluon string bits 
with quarks at the ends or closed rings of gluons, and evaluate matrix elements of the effective 
Hamiltonians between such objects. Since one has a perturbative expression for the operators 
see Eqs. (3.1) and (3.2), one can attempt evaluation of matrix elements between states 
that are constructed in a variety of ways, using quarks and gluons corresponding to different 
scales A. One could even ask if there is a way to calculate a connection between the quark and 
gluon matrices of intermediate widths and reggeized gluon interactions, once one restricts the 
space of states to those that dominate in multi-Regge kinematics. p2 



7. CONCLUSION 

Asymptotically free theories can be analyzed using Hamiltonian approach. The analysis 
can be based on the similarity renormalization group procedure for effective particles. The 
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evaluation of running couplings in the effective Hamiltonians can be carried out without intro- 
duction of wave function renormalization constants and without invoking any properties of the 
vacuum state (in the light-front form of Hamiltonian dynamics). In third order calculations, 
one obtains familiar asymptotic results in scalar 0^ theory, plus an x-dependent finite coun- 
terterm. In QCD, the standard asymptotic freedom form of triple-gluon vertex counterterm is 
supplemented by an ultraviolet diverging and x-dependent counterterm, and by an ultraviolet 
finite a:-dependent counterterm. Effects predicted for QCD by the power counting in k-^ and 

[O] confirmed but the analysis is changed by transition to boost invariant variables k"*" 
and X, and detailed calculations may produce results that are not expected to emerge from 
Feynman diagrams. Mixing between the small x and large k-^ cutoffs indicates a need for a 
new precise definition of the ultraviolet domain in the Hamiltonian approach. Nevertheless, 
one obtains well-defined expressions for effective Hamiltonian interactions without necessity to 
calculate scattering matrix elements for quarks and gluons as if they were observable particles. 

The effective particle calculus preserves cluster properties and allows for evaluation of effec- 
tive Hamiltonians without limitation to any particular set of matrix elements. In other words, 
we can derive integral expressions for matrix elements of effective Hamiltonians in the whole 
Fock space spanned by basis states of effective particles. The renormalization group equations 
are integrated analytically using gaussian similarity form factors and one fully controls off-shell 
behavior of effective vertices that correspond to the initial theory. The effective dynamics is 
invariant with respect to boosts and allows simultaneous analysis of the rest frame and infinite 
momentum frame structure of bound states. 

The effective particle Fock space expansion can converge thanks to the similarity form factors 
in the interaction vertices. The form factors dampen interactions changing invariant masses by 
more than A and thus can tame the spread of eigenstate wave functions for low lying eigenvalues 
into regions of high relative momenta of constituents. This feature may lead to exponential 
convergence of the eigenstate expansion in the effective particle basis. Such convergence is 
not expected in the case of bare particles. The fine structure of effective particles would then 
unfold in the transformation UxM\2 relating effective degrees of freedom at two different scales, 
one corresponding to the binding scale and the other to the high momentum transfer probe in 
question. 

The near-diagonal Hamiltonian matrices allow for transition to new degrees of freedom by 
turning to basis states that are eigenstates of small block matrices on the diagonal. In principle, 
these new degrees of freedom could correspond to mesons and baryons built from constituent 
quarks and gluons, in case of QCD. One can also consider other changes of basis states and 
seek most efficient degrees of freedom, such as strings of gluons with quarks at the ends, for 
solving the nonperturbative eigenvalue problems for narrow effective Hamiltonians. 
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